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We show that topologically protected defect states can exist in open (leaky or lossy) systems even 
when these systems are topologically trivial in the closed limit. The states appear from within the 
continuum, thus in absence of a band gap, and are generated via exceptional points (a spectral tran¬ 
sition that occurs in open wave and quantum systems with a generalized time-reversal symmetry), or 
via a degeneracy induced by charge-conjugation-symmetry (which is related to the pole transition of 
Majorana zero modes). We demonstrate these findings for a leaking passive coupled-resonator opti¬ 
cal waveguide with asymmmetric internal scattering, where the required symmetries (non-hermitian 
versions of time-reversal symmetry, chirality and charge-conjugation) emerge dynamically. 
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Fundamental symmetries appear in a new light when 
they are discussed within the context of open systems, 
where particles escape via leakage to the outside world 
or are absorbed within the material. In these situa¬ 
tions one typically encounters decaying normal modes 
that can be described via the complex eigenfrequen- 
cies uj n of an effective non-hermitian Hamiltonian, with 
H ^ . This description applies, e.g., on the level 

of the Helmholtz equation for dielectric microresonators 
or photonic crystals, where leaky losses enter through 
the boundary conditions while absorption renders the re¬ 
fractive index complex [lj. Such systems still obey reci¬ 
procity, H = H t , while the antiunitary time-reversal 
symmetry THT = H* ^ H T = H is in general broken. 

There has been much recent interest in open settings 
where a generalized antiunitary symmetry VTHVT = H 
still exists EH3- Here V stands for parity, in this con¬ 
text understood to be a unitary involution with V 2 = 1 
that is often realized by a geometric reflection or inver¬ 
sion. True VT symmetry requires amplifying (active) 
regions, arranged such that they are mapped via V onto 
corresponding absorptive parts, while leakage needs to 
be negligible. Eigenfrequencies u> n then are either real 
or occur in a complex-conjugated pair Q, giving rise to 
a novel type of mode competition in lasers US- The 
more general case of a symmetry 


VT{H + i'y'j'PT = H + *7 


(1) 


with a finite offset 7 also encompasses suitably arranged 
passive systems where V now transforms between regions 
of different losses, while uniform leakage is also admit¬ 
ted. The constraints on the spectrum then apply to the 
shifted frequencies f l n = cu n + i'y, which are either real or 
appear along with a partner SI* . Among the many appli¬ 
cations, these features can be used, e.g., to engineer band 
structures in periodic media where the dispersion is still 
effectively real or possesses some well-defined additional 
complex branches TiH- 16 |. 


The advent of topological insulators and superconduc¬ 
tors [ 13 , El has taught us that the classification of uni¬ 
versality classes in terms of time-reversal symmetry is 
incomplete. Two related symmetries, chirality mm 
and charge conjugation [ 21 ], need to be accounted for to 
identify band structures associated with finite topological 
quantum numbers 22 ], with the most prominent conse¬ 
quence being the formation of spatially localized defect 
states at interfaces (points, edges or surfaces) between 
topologically distinct domains. In particular, a unitary 
chiral anti-symmetry X(H — = — ( H — fl (°)) en¬ 

forces the spectrum to be symmetric around a central 
frequency giving rise to frequency pairs ± f l n . 
A similar spectral constraint is also enforced by the an¬ 
tiunitary charge-conjugation symmetry C = TX that 
can appear in superconducting systems; this yields pairs 
fl(°) + — fi* and stabilizes unpaired resonances 

(broadened Majorana zero modes) at Refl„ = 0 [26l — 
1281 ]. Within this conventional classification, a necessary 
requirement for topological nontriviality is the existence 
of a band gap, into which the defect states then fall. 
The nascent field of topological photonics 2^- 441 has 
embarked to realize photonic analogues of these sym¬ 
metries, while the (beneficial or detrimental?) role of 
non-hermitian loss and gain has only been considered in 
settings which are already topological in the hermitian 
limit MM- 


Here we identify, for the simple example of a coupled- 
resonator chain, a mechanism by which topologically pro¬ 
tected defect states can appear in open (non-hermitian) 
systems even when their closed (hermitian) limit is topo¬ 
logically trivial. The defect states form at an inter¬ 
face of two regions with different non-hermiticity, and 
appear from the continuum of the band structure via 
two distinct symmetry-protected spectral transitions— 
where one is associated with PT-symmetry (thus, non- 
hermitian time-reversal symmetry), while the other is as¬ 
sociated with PC-synnnetry (thus, the analogously gener- 
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FIG. 1. (color online) (a) Coupled-resonator waveguide with 
internal asymmetric scattering (couplings A and B) between 
a counterclockwise (CCW) and a clockwise (CW) wave com¬ 
ponent, and coupling W between CCW and CW components 
in neighboring resonators. For real couplings this system re¬ 
alizes non-hermitian versions of time-reversal symmetry, chi¬ 
rality and charge conjugation, (b) Phase diagram for the bulk 
dispersion ©. The dispersion can be real and gapped, exhibit 
2 or 4 exceptional points (EPs) at which real and imaginary 
branches meet, or be fully imaginary. In the examples in (c), 
A/W ~ 4 and B/W = 3 (real gapped dispersion), A/W = 4 
and B/W = 1 (2EPs), A/W = 1 and B/W = -1 (4EPs), as 
well as A/W = 4 and B/W = —3 (imaginary dispersion), (d) 
Corresponding representation of the complex bands on the 
Bloch sphere (upper band Q+, red; lower band S2_, blue). 


alized charge-conjugation symmetry that arises from the 
simultaneous presence of a chiral symmetry). Therefore, 
robust states as desired for topological photonics can be 
obtained by combining symmetries with non-hermitian 
effects that go beyond the universal classification of elec¬ 
tronic systems. 

Dynamical realization of non-hermitian parity-time 
and charge-conjugation symmetry .— We first describe 
how the required symmetries can be implemented in a 
simple passive resonator chain, with losses solely pro¬ 
vided due to leakage but without any need of absorp¬ 
tion or amplification (which then also translates to anal¬ 
ogous open quantum systems). This can be achieved 
in a coupled-resonator optical waveguide (CROW) |52f 
Ib3 | which consists of identical asymmetric cavities [55|, 
as sketched in Fig. Oja). Each individual resonator fea¬ 
tures two modes—a counterclockwise (CCW) propagat¬ 
ing mode with amplitude a n and a clockwise (CW) prop¬ 
agating mode with amplitude b n , which we group into a 
vector ifn = ( a n ,b n ) T . Resonator arrays in which these 
two modes are well decoupled feature in setups that re¬ 
alize photonic topological edge states in analogy to the 
quantum-Hall effect 34| 381. In our setting, however, the 


internal coupling between these modes is desired, and the 
key feature is that this coupling can be made asymmet¬ 
ric by opening the system, even if no magnetic field is 
applied—as has been established in recent works on indi¬ 
vidual resonators 55-6Cj. The coupling of the modes is 
then described by a nonhermitian internal Hamiltonian 

mm 


h = 


A 

B 


( 2 ) 


where the constants A ^ B* and account for the 
asymmetric internal scattering and the losses within the 
cavity. Throughout the chain, the coupling between ad¬ 
jacent resonators is dominantly between CCW and CW 
waves, so that the coupling matrix is SE3 


t = 


0 W 
W 0 


( 3 ) 


In coupled-mode approximation, the stationary wave 
equation then takes the form 


U)4> n = hlfn + t(lp n+ 1 + 4>n-l), 


( 4 ) 


which admits Bloch solutions %l> n = exp (ikn)^. The as¬ 
sociated Bloch Hamiltonian is 

i \ ( n (0) A + 2W cos(k) \ 

m =(B + 2Wcosm «<»> )’ (5) 

and leads to the dispersion relation u>±(k) = + 

Ofc(fc), 


D±(k) = ±yJ{A + 2W cos k)(B + 2 W cos k ), ( 6 ) 

where the subscript ± labels two bands. The symmetry 
about is a consequence of a chiral symmetry with 
Xif n = a z if n , which maps the bands onto each other. 
The chiral symmetry is thus realized by the freedom of 
the relative sign of the CCW and CW amplitudes (a 
gauge freedom compatible with time-reversal symmetry 
which generically appears in systems with two mutually 
coupled sublattices). 

As shown by exact numerical calculations [60, |(jll , for 
representative resonator geometries A, B and W are al¬ 
most real, and can be further tuned towards real values 
by adjustments of a few shape parameters. We thus ne¬ 
glect the imaginary parts of these parameters. Apart 
from an offset 7 = ilmfl( 0 ), the dispersion is then ei¬ 
ther real (in some range of k) or purely imaginary (in the 
complementary range of k). These k ranges are joined by 
degeneracies, known as exceptional points 62|j66j. where 
w±(fc) = thus cos k = —A/2W or cos A; = —B/2W. 
A completely real dispersion with a gap is achieved if 
|A/2Wj > 1 and \B/2W\ > 1, provided AB > 0 [see the 
phase diagram in Fig. [Ub) and representative dispersions 
in Fig. [He)]. 
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The underlying symmetry can be made explicit by a 
basis change, </> n = 2~ 1 / 2 (i<r x + t)if n where <j x is a Pauli 
matrix, after which the Bloch Hamiltonian takes the form 

r,,, _ 1 / 2Q^ + i(A-B) A + B + 4Wcosk\ 

h ^ } ~ 2 ^ A + B + 4Wcosk 2 QW-i(A-B) )' 

(7) 

The Hamiltonian is now symmetric, as required by reci¬ 
procity (which is hidden in basis of CCW and CW modes 
since propagating waves are complex), and furthermore 
exhibits a passive VT symmetry (jXJ) with V = a x , and 
7 = —Imfl*°l The chiral symmetry is transformed to 
X = a y , and commutes with VT. This realizes all 
the symmetries mentioned in the introduction, including 
VC = VTX , with respect to the central frequency Cl^ . 
From here on, we work in terms of the shifted frequencies 
Cl = to — for which the dispersion is directly given 
by Eq. d6j). 

Having established these symmetries we now return to 
the basis of CCW and CW modes and discuss topolog¬ 
ical aspects of the band structure. For this we consider 
the k dependence of the Bloch vectors 'F(fc), which we 
interpret as pseudospins with polarization vector P = 
((a x ,a y ,a z )). In the hermitian limit B = A (both real), 
the two bands Cl±(k) = ±(A + 2W cosk) arise from 
k -independent pseudospins 'k-t = 2 _1 / 2 (1, ±1) T , with 
P = (±1,0,0) pointing along the x axis. The absence 
of any winding of the pseudospin renders the system 
topologically trivial, so that we do not expect any de¬ 
fect states in the presence of interfaces, even if there is a 
gap. In the non-hermitian case, we can write ’P±(fc) oc 
(A+2W cosfc, Cl±{k)) T . As shown in Fig.OJd), the polar¬ 
ization vector now acquires k dependence; it is confined 
to the xz plane when the dispersion is real and to the yz 
plane when the dispersion is imaginary. These branches 
are again joined at the exceptional points, where P points 
up or down along the z axis, with P-(k) = R z (Tr)P + (k) 
related by a 7r rotation about the z axis. In particular, 
the way these points are connected depends on whether 
A > B or A < B (with the two cases related by a rota¬ 
tion R x (tt) by 7T about the x axis). Does the system now 
admit defect states? 

Defect states .— In order to answer this question, we 
create a defect in the chain by inverting the orientation of 
the resonators in half of the system [see Fig. [51(a)]. From 
the traditional perspective of hermitian systems, the de¬ 
fect cannot be classified as topological, and does not give 
rise to any defect states. In the non-hermitian setting, we 
will see that the defect acquires topological features in a 
spectral phase transition at which localized defect states 
emerge [as illustrated in Fig. Efb)]. The phase transition 
takes the form of a PT-induced exceptional point along 
one part of the phase boundary, while it is associated 
with a PC-induced degeneracy along the other parts of 
the phase boundary [this is summarized in Figs. [5][c,d), 
to which we refer throughout the remaining discussion]. 
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FIG. 2. (color online) (a) Coupled-resonator waveguide with 
a defect, created by inverting the orientation of the resonators 
in half of the system. In the closed limit, the system is trivial, 
and the defect does not create any bound states, (b) Defect 
states in a system of 300 resonators, with A/W = 1.9 (left 
panel) and A/W = 2.5 (right panel), while B/W = 1. (c) 
Phase diagram indicating the existence of defect states, as 
well as their extended-state precursors (realizing perfect in¬ 
terband transitions, PIC). The boundaries of the defect phase 
are given by degeneracy conditions. At the PT boundary, 
the extended PIC states bifurcate into pairs Cl„, of de¬ 
fect states that are related by the VT symmetry. At the PC 
boundary one encounters a degeneracy of charge-conjugated 
partner states — fiji, beyond which the defect states are 

non-normalizable. (d) Bloch-sphere position of the defect 
states (rods) relative to the bulk dispersion (lines) [parameters 
A/W = 1.7, 1.81 (PT), 1.9, 2.5, 3.32 (PC), 4 with B/W = 1, 
as indicated by the white circles in (c)]. 
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It is easy enough to identify the conditions for the for¬ 
mation of defect states. In the presence of the defect, the 
wave equation takes the form 

ttlpn = Klpn + t(^} n+ 1 + Pn- 1), (8) 

where now h n = (h — for n < 0 (left half of the 

chain) and h n = {h — fV°)) T for n > 0 (right half of the 
chain). At any fixed f2, in each half of the system the 
solutions are still obtained from a superposition of Bloch 
waves with 


2 W cos k± = c± = 


A + B 


± 


(A-B) 2 


+ ^ 2 , (9) 


associated with pseudospins 


*±’(0) o< ( ^ q C± ) . *?(«) « ( ^ ) ■ 

( 10 ) 

Each of the values c+, c_ are associated with a pair of 
Bloch waves with propagation factors exp(ifc±) = A±, 
exp(— ik±) = (A±) _1 , where we choose k± such that 
|A±| > 1 if k± is complex. We now match the solutions 
with propagation factor exp(ifc±) in the left part to the 
solutions with propagation factor exp(— ik±) in the right 
part. One then finds the condition 


(n-2W) 2 n-ABn = (A-b) 2 w/2 ( 11 ) 


for defect states with a symmetric wavefunction, and 


states. The same scenario occurs simultaneously for the 
chirality-related solutions of Eq. 03 . The second kind 
of transition appears along the curves labeled PC, where 
A 2 + GAB + B 2 = 32W 2 . There, a complex solution il n 
of Eq. (fill) coalesces with a charge-conjugated solution 
ofEq. H3, meaning that they are purely imaginary. 
This is similar to the pole transition of broadened Majo- 
rana zero-modes, which are then pinned to the imaginary 
axis and become their own charge-conjugated partner. 
These transitions also occur in skew-Hamiltonian ensem¬ 
bles governing the topological transitions in Josephson 
junctions 67J. In the present problem, the PC transition 
signals the point where the matching conditions can only 
be fulfilled by combining decaying with increasing wave- 
functions, which occurs when one of the wave numbers 
k .|_, k- crosses the real axis. 

At both types of transition, the defect states there¬ 
fore interact with the real branch of the dispersion re¬ 
lation f l±(k) (for the exceptional points along the PT 
boundary), or with the purely imaginary branch of this 
dispersion relation (for the charge-conjugation-induced 
degeneracy along the PC boundary). On the level of 
the wavefunctions, this interaction is again revealed via 
the corresponding polarization vectors. Focussing on the 
wave function in the left part of the system, we have 


p[ L \n) = -Pl L \n*) ( 13 ) 

for the two partial waves given in Eq. m, while 
the VT and chiral symmetries relate 


- (n + 2W) 2 n + ABn = (A - b) 2 w/2 ( 12 ) 

for defect states with an antisymmetric wavefunction. 
Because of the VT symmetry, solutions are again either 
real or appear in complex conjugated pairs. Further¬ 
more, chirality maps fl —>- —f2, corresponding to a trans¬ 
formation between symmetric and antisymmetric wave- 
functions. Thus, the defect states come in a quadruple of 
frequencies f l n , fi*, —fi n , —fi*. The condition that the 
corresponding wave function indeed decays leads to the 
phase diagram in Fig. He). 

In this phase diagram, the hermitian case A = B 
defines a diagonal. The defect states are confined to 
a region away from the diagonal, which is bounded by 
two different transitions. Along the curves labeled PT, 
where 27(A + B ) 4 = 16A 2 B 2 (1 + AB/W 2 ) + 8(8 W 2 + 
9AB)(A + B) 2 , a pair of real solutions fl n of Eq. (fill) 
bifurcates into a pair of complex-conjugated solutions. 
Before this exceptional point, the solutions are real, with 
|A±| = 1, and describe the scattering of an incoming ex¬ 
tended state in one band into an outgoing extended state 
in the other band. This region of perfect interband con¬ 
version is labeled PIC. At the exceptional point, the prop¬ 
agation factors A.), of the two solutions coalesce, and so 
do the factors A_; beyond the exceptional point we then 
have |A±| > 1, giving rise to properly normalizable defect 


Pi L \n*) =-R v (n)Pi L \n), (14) 

Pi L \-n) = R z (n)Pi L \n), (15) 

where R a denotes a n rotation about axis a. [We also 
have P^ i \fl) = R x (n)Pi L \n).] In Fig. Hd), we show 
how these polarization vectors interact. In the PIC 
phase, each state corresponds to a pair of opposite vec¬ 
tors P^ confined to the xz plane (the locus of the real 
dispersion branch), with chirality-related partner states 
connected by a n rotation about the z axis. At the PT 
transition, the vectors bifurcate and move out of the xz 
plane. In the defect phase, the two vectors P^ for a 
given defect state make an angle, but remain related by 
a 7r rotation about the y axis. The partner state with fre¬ 
quency kl* points into the opposite direction, while the 
chirality-related states are still obtained by a 7r rotation 
about the z axis. At the PC transition, each state col¬ 
lides with a charge-conjugated partner at a point in the 
yz plane (the locus of the imaginary dispersion branch). 
These interactions all occur at symmetry-protected posi¬ 
tions, which renders the defect phase topologically stable. 
For numerical verification of this robustness in finite and 
disordered systems see [§§]. 

In summary, robust defect states can exist in open sys¬ 
tems that are topologically trivial in the closed limit. 
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We illustrated this for a leaky optical resonator chain 
where defect states appear at an interface between re¬ 
gions in which hermiticity is broken in different ways 
(in contrast, non-hermiticity is not sufficient to create 
edge states at the end of a finite sample). The states are 
topologically protected as they arise in spectral phase 
transitions that are linked to the spontaneous breaking 
of fundamental symmetries (parity-time symmetry and 
charge-conjugation symmetry) for a sufficient degree of 
non-hermiticity in the system. The required symmetries 
are realized when the couplings in the propagating wave 
basis are real. As this does not require any absorption, 
the formation mechanisms described here also translate 
to analogous geometrically open quantum systems, in¬ 
cluding electronic systems which are suitably coupled to 
external reservoirs. Our observations raise new ques¬ 
tions, such as whether it is possible to characterize these 
systems in terms of topological quantum numbers, and 
more generally whether they can be understood by a suit¬ 
able extension of the conventional topological classifica¬ 
tion of closed systems. 

We thank Jan Wiersig for useful discussions and com¬ 
ments. This research was supported by EPSRC via grant 
EP/J019585/1. The data created during this research is 
openly available (7lj| . 
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Appendix: Finite systems and disorder 

To illustrate the robustness of the defect states and the 
PT and PC phase transitions, we extend our considera¬ 
tions to finite systems as well as situations in which A 
and B are subject to further asymmetries and disorder, 
with different fixed or average values in the left and right 
halves of the system. The left-right symmetry (which we 
did not exploit in our general considerations) is then bro¬ 
ken, while the VT and VC symmetries continue to hold in 
any given system. To facilitate this discussion, we first 
introduce a number of useful quantities to characterize 
and identify the defect states in such general settings. 


Figure lAlf a) shows a defect state in a finite system 
of 300 resonators, with A — W and B = 2 W. The 
defect is placed into the center, so that there are two half 
chains of 150 resonators. The amplitude a n refers to the 
CCW component in the left half of the chain and the CW 
component in the right half of the chain, while b n refers 
to the CW component in the left half of the chain and 
the CCW component in the right half of the chain; see 
coloring of sites in Fig. 2(a) in the main text. 

By changing the couplings smoothly according to B = 
x , we can steer the system from the defect phase over 
the PT and PC phase boundaries. The trajectories of 
the quadruple of defect eigenvalues and their precursors 
are shown in Fig. lAlf b). At the PT transition, we find 
that the degeneracy on the real axis is exactly realized 
even in the finite system. We relate this to the fact that in 
the wave matching problem, the defect state then morphs 
into the PIC precursor state, which simply is a specific 
example of an extended state. At the PC transition, in 
the wave-matching picture the defect state morphs into 
a non-normalizable state. As the spectrum of the finite 
system cannot lose an eigenvalue, the exact degeneracy 
on the imaginary axis is then lifted; we find that the 
state hybridizes with extended states as soon as its decay 
length becomes comparable with the system size. The ex¬ 
act degeneracy is then approached more and more closely 
when the system size is increased. 

In Fig. lAlf c). we quantify the closeness to spectral 
degeneracy by the Petermann factor 


(L\L)(R\R) 

\(L\RW 


(A.l) 


This is a measure of non-orthogonality of the left eigen¬ 
states (L | and right eigenstates | R) which diverges at 
non-hermitian degeneracies [65l . l66| (note that the no¬ 
tion of left and right eigenstates does not refer to the left 
and right half of the chain, but to the left and right mul¬ 
tiplication of the states with the non-hermitian effective 
Hamiltonian). The Petermann factor also quantifies the 
sensitivity of the system to general perturbations [66}, in¬ 
cluding quantum noise 69|,|70|, which thus can be used to 
map out the boundaries of the defect phase in such gen¬ 
eral settings. We find that the Petermann factor becomes 
large at the PT and PC phase boundaries, in accordance 
with the spectral degeneracies that one approaches there. 

In order to assess the localization of the defect state 
we use the inverse participation ratio 


IPR = ^(|a„| 4 + |6„| 4 ), (A.2) 

n 

which is shown in Fig. lAlI d). As expected, the IPR is 
large within the defect phase. 

As shown in Fig. Ella) for an asymmetric system 
(A = W and B = 2W in the left half and A = 0.5 W, 
B = 2.5 W in the right half of the chain), well-defined de¬ 
fect states still exist in this more general situation, where 
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FIG. Al. (color online) (a) Mode profile of a defect state in 
a finite system consisting of 300 asymmetric resonators with 
A = W and B = 2W [blue: intensity |a„| 2 of the CCW com¬ 
ponent in the left half and CW component in the right half; 
red: intensity |6 ra | 2 of the CW component in the left half and 
CCW component in the right half]. Panel (b) shows the tra¬ 
jectories of the eigenvalues associated with the defect states 
and their predecessors as the value of B = x is changed over 
the range 1 < x < 4.5 at fixed A = W. (c) Petermann factor 
K, Eq. ED, as a function of x. The Petermann factor be¬ 
comes large at the PT and PC phase boundaries, as well as at 
the exceptional points in the PIC phase, in accordance with 
the spectral degeneracies encountered there, (d) Inverse par¬ 
ticipation ratio (IPR), Eq. (1X2)1 . showing strong localization 
of the defect state in the expected range. 


we now observe a different decay of the intensity into the 
two half chains. By changing the couplings smoothly, 
now according to B = x (left half) and B = x + 0.5 W 
(right half), we again find that the defect phase termi¬ 
nates at PT and PC transitions, where the Petermann 
factor becomes large. The IPR is again large within the 
defect phase. 


Figure [A3] displays the same information in the pres¬ 
ence of additional disorder, where the couplings A and B 
in each resonator are modified by independent perturba¬ 
tions y from a uniform distribution y £ [—0.1 W, 0.1 W\. 
This breaks the translational invariance in each half chain 
an induces Anderson localization, with a localization 
length that competes with the defect state, especially 
when the boundaries of the defect phase are approached. 
While the disorder shifts the PT and PC transitions in 
parameter space, the defect phase remains well defined, 
and the features of the defect states are remarkably ro¬ 
bust. 













FIG. A2. (color online) Same as Fig. IA1I but for couplings 
A = W and B = 2 W in the left half and A = 0.5 W, B = 
2.5 W in the right half of the chain (panel a), as well as A = 
W, B = x (left half) and A = 0.5 W, B = x + 0.5 W (right 
half) in the remaining panels. 



Resonator number n 




FIG. A3, (color online) Same as Fig. IA2I but with the cou¬ 
plings A and B in each individual resonator perturbed by 
independent amounts y, drawn from a uniform box distribu¬ 
tion y G [—0.1 W, 0.1 W]. The results represent one disorder 
realization in panel (a), while the other panels compare the 
results for three disorder realizations. 





















